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ABSTRACT 
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We report the detection of a planetary system around BD-08°2823, that includes at least one Uranus-mass planet and one Saturn- 
mass planet. This discovery serendipitously originates from a search for planetary transits in the Hipparcos photometry database. 
This program preferentially selected active stars and did not allow the detection of new transiting planets. It allowed however the 
identification of the K3V star BD -08°2823 as a target harboring a multiplanet system, that we secured and characterized thanks to an 
intensive monitoring with the HARPS spectrograph at the 3.6-m ESO telescope in La Silla. The stellar activity level of BD -08°2823 
complicates the analysis but does not prohibit the detection of two planets around this star. BD -08°2823b has a minimum mass of 
14.4 ±2.1 M e and an orbital period of 5.60 days, whereas BD -08°2823c has a minimum mass of 0.33 ± 0.03 Mj up and an orbital 
period of 237.6 days. This new system strengthens the fact that low-mass planets are preferentially found in multiplanetary systems, 
but not around high-metallicity stars as this is the case for massive planets. It also supports the belief that active stars should not be 
neglected in exoplanet searches, even when searching for low-mass planets. 

■ Key words, planetary systems - techniques: radial velocities - stars: individual: BD -08°2823 
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Q\ I 1. Introduction hitherto plays a major role in the improvement of the knowledges 

. of exoplanets. 

•' The HARPS spectrograph (Mayor et al. I2003I I is operating since We report here the detection of a new multiplanetary system. 
.£h ! 2003 at the 3.6-m ESO telescope in La Silla, Chile. This is It serendipitously originates from the radial velocity follow-up 
^ . a fiber-fed, environmentally stabilized high-resolution echelle accompanying a search for planetary tra nsits i n the Hipparcos 
$_i ; spectrograph dedicated to high-precision radial velocity mea- epoch photometry annex (Perryman et al. fTWfl . Despite its sig- 
\ surements. Thanks to the exoplanetology programs that are con- nificant stellar activity, the star BD -08°2823 was identified dur- 
ducted with it, it allowed numerous extra-solar planets studies i ng this follow-up as a possible host of extra-solar planets. An 
and discoveries, in the r anges of Jupiter-mass planets (e.g. Pepe intensive monitoring of its radial velocities was subsequently 
et al. 2004 Moutou et al.|2009a) , low-mass planets (e.g. Lovis et carried out with HARPS, that led to the detection of a plane- 
al.|2006 Mayor et al. 2009), planets around early- (e.g. Desort tary system around this star, including at least two planets. We 
et al. I2008I Lagrange et al. [2009 ) or late-type stars (e.g. Bonfils briefly present the search for planetary transits in the Hipparcos 
et al. I2007I F orveil le et al. I2009), or t ransiting planets (e.g. photometry in § The HARPS observations of BD -08°2823 
Bouchy et al.. I2008I Queloz et al. |2009l >. The essential quality are presented in § [3] together with the stellar properties of this 
of HARPS (the High Accuracy Radial velocity Planet Searcher) target. The planetary system is characterized in § @J then the re- 
is its high stability, which results in a sub-ms -1 accuracy in the suits are discussed in §[5] 
radial velocity measurements, on time-scales of several years. 
This allows the detection of planets in the Neptune or Super- 
Earth mass ranges, on progressively increasing orbital periods as 2. Search for transiting planets in the Hipparcos 
the time-span of the monitoring is growing up. This instrument database 

The goal of our initial search in the Hipparcos database was to 

* Based on observations made with HARPS spectrograph on the 3.6- tr y to find new transiting planets, especially around bright stars, 

m ESO telescope at La Silla Observatory, Chile, under the programs ID Transiting planets could allow numerous important studies to 

072.C-0488, 074.C-0364 and 078.C-0044. be performed. These studies include planetary radii, masses and 
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densities measurements, accurate determination of the inclina- 
tions and the eccentricities of the orbits, possible detections of 
transit timing variations due to additional companions, measure- 
ments of the spin-orbit (mis)alignment angles, detections of the 
emission or absorption planetary spectra, or even satellites or 
rings detections. The famous planets transiting in front of the 
bright stars HD 209458 and HD 189733 are those that allow the 
most accurate measurements and the deepest investigations. 

The Hipparcos epoch photometry annex contains between 
~ 40 and ~ 300 photometric measurements performed during 
the 1226-day duration of the mission for each of the 118 204 
stars of the catalog, in the magnitude limit V < 10. With about 

0. 1% chance that a given star harbors a transiting extra-solar 
planet, the Hipparcos all-sky survey must contain photomet- 
ric measurements for tens of transiting hot Jupiters. This be- 
lief is reinforced by the two cases of transiting planets a pos- 
teriori re-discovered in Hipparcos data: HD 209458b (Robichon 
& Arenou|2000l Castellano et al. [20001 Soderhielm [19991 and 
HD 18 9733b (Bouchy et al. [20031 Hebrard & Lecavelier des 
Etangs 2006). Thanks to the long time baseline, these data also 
allow the orbital period to be measured with an accuracy of the 
order of the second. 

Up to now, there are the only two transiting planets that have 
been found in the Hipparcos database. Indeed, with photometric 
variations of about 1 % or less, such transits are difficult to iden- 
tify in the Hipparcos data, that present a mean photometric accu- 
racy of the same order of magnitude. The poor time coverage, by 
comparison with dedicated surveys as SuperWASP, XO, CoRoT 
or Kepler, is another difficulty. Jenkins et al. ([2002 ) concluded 
that because of its poor photometric quality, the Hipparcos cata- 
log does not represent a likely place to detect planets in the ab- 
sence of other informations. It might however provide planetary 
transit candidates for follow-up observations. Laughlin (2000) 
has searched in the Hipparcos epoch photometry for transit- 
ing planets within 206 metal-rich stars. None have been con- 
firmed thereafter. 

We present in Appendix [A] the systematic search we man- 
aged in the Hipparcos epoch photometry annex for periods com- 
patible with planetary transits. We constructed a ranked list of 
candidates for follow-up radial velocity measurements, based on 
the depth and significance of the planetary transits that could 
be detected in the Hipparcos photometry. The simulations we 
performed indicated that our detection rate for transiting planets 
is low, of the order of 2 %; so the chances to find new transit- 
ing planets in the Hipparcos data are small (see Appendix [At. 
However, the particularly high interest for the potential discov- 
ery of a planet transiting a bright star pushed us to perform a 
radial velocity follow-up of our candidates. 

We observed 194 of these selected, ranked targets with 
HARPS, in order to search for radial velocity variations in agree- 
ment with the transiting candidates found in the Hipparcos pho- 
tometry. The HARPS setup, the spectra extraction and the ra- 
dial velocity measurements were identical to those described 
below in § 13.11 Most of these observations were performed in 
December 2004 as part of the program 074.C-0364. Radial ve- 
locity variations larger than 20 m s _1 were measured for 37 stars, 

1. e. 19% of our observed sample. Those variations are mainly 
caused by stellar activity. Indeed, most of the targets with large 
radial velocity variations show emissions in the Ca n H & K lines 
(3934.8 A and 3969.6 A). The level of these emissions is quanti- 
fied with the activity S-index (Mount Wilson System), which is 
converted to the logR' HK index (Santos et al. 2000). The majority 
of the variable targets have log R' RK indexes larger than -4.7, in- 



dicating prominent chromospheric activity; this can explain such 
level of radial velocity variations (Santos et al.[2p00). 

Some targets first seemed to exhibit periodic radial velocity 
variations, but subsequent monitoring with HARPS showed that 
these variations were not persistent with time, in agreement with 
transient activity processes on the stellar surface due to activity 
(flares, spots, plages...) modulated by the stellar rotation. In ad- 
dition, in most of the cases, analysis of the line profiles using 
bisectors of the cross-correlation functions shown variations in 
the shape of the lines with time. This indicates that the observed 
radial velocity variations are not due to Doppler shifts of the 
lines, but rather to deformations of the shape of the spectral lines. 
Some cases revealed a clear anti-correlation between radial ve- 
locity and line-bisector orientation, which could be understood 
as the signature of cool spots on the stellar photospheres (see, 
e.g., Queloz et al I200T1 Melo et al. [20071 Desort et al. [20071 
Boisse et al. |2009b . 

Thus, our procedure seems to preferentially select active 
stars. This should be contrasted with transiting candidates ob- 
tained from the photometric surveys dedicated to transiting 
planet searches, which are mainly eclipsing binaries or transit- 
ing planets (see, e.g., Pont et al. 2005 ). None of the targets that 
we observed through this program have shown radial velocity 
variations in agreement with a reflex motion due to a hot Jupiter. 
Instead, most of the radial velocity variations seem to be caused 
by stellar activity. 

However, these observations allowed the serendipitously dis- 
covery of a new planetary system, without transit detection in 
the Hipparcos data. This system orbits BD -08°2823, which 
is an active star, but among those with the lowest logR' UK in- 
dexes within our monitored targets. This activity level compli- 
cates planet detection but does not prohibit it in this case. We 
concentrate below on this target and the detection of its harbored 
planetary system. 

3. Observations and properties of BD -08°2823 

3.1. HARPS observations of BD -08°2823 

We observed BD -08°2823 using the HARPS spectrograph at 
the 3.6-m ESO telescope in La Silla. The bandpass of the spec- 
tra ranges from 3800 A to 6900 A, and the resolution power is 
R = 115 000, with a fiber diameter of 1 arcsec. The spectra were 
extracted from the detector images with the HARPS pipeline, 
that includes localization of the orders on the 2D-images, opti- 
mal order extraction, cosmic-ray rejection, wavelength calibra- 
tion and corrections of flat-field. The pipeline then performs a 
cross-correlation of the extracted spectra with a numerical mask 
(K5-type in this case), and finally measures the radial velocities 
from Gaussian fits of the cross-correlation functions (CCFs), fol- 
lowing the method described by Baranne et al. ( 1996) and Pepe 
et al. (120021) . The full dataset we use for BD -08°2823 includes 
83 spectra. All the exposures (but the first one) were obtained 
without simultaneous thorium calibration. The exposure times 
for BD -08°2823 range between 4 and 15 minutes, allowing a 
signal-to-noise ratio per pixel between 30 and 80 to be reached 
around 550 nm. The total exposure time is about 1 1 hours. 

The HARPS radial velocities of BD -08°2823 as a func- 
tion of time are plotted in the upper panel of Fig.Q] They span 
1841 days, corresponding to 5.0 years. A few measurements 
were obtained in late-2004 as part of the program 074.C-0364 
(see § [2]), then the target was regularly monitored. A few extra 
measurements were acquired in late-2006/early-2007 as part of 
the program 078.C-0044, but most of the data were secured as 
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part of the Guaranteed Time Observations (GTO) survey pro- 
gram 072.C-0488 (Mayor et al. 120031 . This is the case in par- 
ticular for the intensive series in 2007-2008 (BDJ-2 400 000 be- 
tween 54150 and 54650). There was indeed an agreement to 
follow, as part of the GTO time, a few promising targets identi- 
fied from the Hipparcos original program. 
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Fig. 1. Top: Radial velocity HARPS measurements of 
BD -08°2823 as a function of time (see § |3), and Keplerian fit 
with two planets (see §|4]l. The orbital parameters corresponding 
to this fit are reported in Table [3] Bottom: Residuals of the fit 
with 1-cr error bars. 



Table 1. HARPS measurements of BD -08°2823 (full table 
available electronically)^. 



BJD 


RV 


cr(RV) 


FWHM 


Bis. 




-2400 000 


km s" 1 


km s -1 


km s 


km s" 1 




53007.83343 


53.3353 


0.0020 


6.281 


0.0222 


-4.702 


53358.83576 


53.3408 


0.0016 


6.248 


0.0183 


-4.731 


53359.80161 


53.3479 


0.0016 


6.255 


0.0154 


-4.725 


53360.81961 


53.3563 


0.0015 


6.252 


0.0226 


-4.742 


54648.45269 


53.3538 


0.0012 


6.322 


0.0326 


-4.709 


54813.83365 


53.3511 


0.0015 


6.291 


0.0273 


-4.755 


54840.78987 


53.3562 


0.0010 


6.300 


0.0294 


-4.732 


54848.79967 


53.3565 


0.0010 


6.283 


0.0257 


-4.751 



f: the six columns are respectively the Bary centric Julian Day of the 
observations, the radial velocity and its uncertainty, the full width 
at half maximum of the CCF, the span of its bisector, and the Can 
activity index. 



(n = 23.76 ± 1.61 mas) implies a distance of 42.2 + 2.9 pc. The 
Hipparcos color is B - V = 1.071 + 0.010. From spectral anal- 
ysis of the HARPS data using the method presented in Santos et 
al. (|2004al l, we derive the temperature T eff = 4746 + 63 K, the 
gravity logg = 4.13 + 0.26, the metallicity [Fe/H] = -0.07 + 
0.03, and an uncertain age of 4.5 + 4.0 Gyr. The stellar mass 
we obtain is M* = 0.74 M , with a formal error bar of +0.2 M . 
Following Fernandes & Santos ( 120041 1. we rather adopt a conser- 
vative + 10 % uncertainty, corresponding to +0.07 M . We derive 
a projected rotational velocity v sin z* = 1 .4 km s from the pa- 
rameters of the CCF using a calibration similar to that presented 
by Santos et al. (120021 1. Acoustic oscillations are not averaged 
out in the data, as there are on timescales similar or shorter than 
the exposure times we used. However, for a K-type star, oscil- 
lation amplitudes are expected to be negligible by comparison 
with stellar activity effects on the measured radial velocities. 



The derived radial velocities are reported in Table [T] The 
accuracies are between 1.0 and 2.4 ms _1 , typically around 
1.6 ms~'. This includes photon noise but not the jitter due to 
stellar oscillations or activity (§ |3.2| >. Wavelength calibration and 
spectrograph drift uncertainties are negligible with respect to the 
photon noise. Four exposures with radial velocity uncertainties 
larger than 2.4 ms _1 were not included in the final dataset of 
83 spectra. 

The radial velocities show a significant ll-ms -1 disper- 
sion (43 ms _I peak-to-peak), well over the expected accuracy 
(Fig. [1] upper panel). The CCF from which those radial veloci- 
ties were measured show parameters that also significantly vary 
with time. Their full widths at half maximum fluctuates between 
6.22 + 0.01 and 6.33 + 0.01 kms -1 , and their contrasts from 
46.80 + 0.05 and 47.50 + 0.05 % of the continuum. The bisector 
span of the CCF also show a dispersion, at the level of 7 ms _1 
(see Fig. [2] upper panel). Thus, Doppler shifts of the spectral 
lines of BD -08°2823 are not necessarily the only explanation 
for the observed radial velocity variations; changes in the shape 
of the lines, as those due to stellar activity, are likely to be the 
cause of at least a part of those variations. 

3.2. Stellar properties and activity 

Table [2] summarizes the stellar parameters of BD -08°2823 
(HIP 49067, SAO 137286). According to the SIMBAD database, 
this is K3V star of magnitude V = 9.86. Its Hipparcos parallax 



Table 2. Stellar parameters for BD -08°2823. 



Parameters 


Values 


m v 


9.86 


Spectral type 


K3V 


B - V 


1.071+0.010 


Parallax [mas] 


23.76+ 1.61 


Distance [pc] 


42.2 ± 2.9 


vsini* [kms -1 ] 


1.4 


P ml [days] 


26.6+1.5 


logtf HK 


-4.8 + 0.1 


[Fe/H] 


-0.07 ± 0.03 


Teff [K] 


4746 ± 63 


l0g£ [Cgs] 


4.13+0.26 


Mass [M G ] 


0.74 ± 0.07 


Age [Gyr] 


4.5 + 4.0 



The core of the large H & K Can absorption lines of 
BD -08°2823 show emissions, indicating chromospheric activ- 
ity. We adopted the averaged value logR' HK - -4.8 ±0.1, but 
it significantly varies with time, between extrema logR' HK = 
-4.68 ± 0.01 and -4.88 ± 0.01. Such stellar activity would im- 
ply a significant jitter on the observed stellar radial velocities. 
For a K-type star with this level of activity, Santos et al. (2000) 
predict a dispersion up to 10 ms _1 for the stellar jitter. This is 
the order of magnitude of the dispersion of our measurements 
of BD -08°2823 (§EB. A classic method to identify stellar ac- 
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tivity as the main cause of radial velocity variations is to look 
for anti-correlation between radial velocity and the bisector span 
(see, e.g., Queloz et al. 2001 Boisse et al. 120091 . The upper 
panel of Fig. [2] shows the bisector spans as a function of the 
radial velocities for BD -08°2823. As seen above, the bisector 
is significantly varying; however, no clear relations are seen be- 
tween the two signals. Yet, this technic is less sensitive for stars 
with slow rotations, as this is the case for BD -08°2823. Indeed, 
according to Noyes et al. (1984) and Mamajek & Hillenbrand 
( 120081 , the activity level for BD -08°2823 implies a stellar rota- 
tion period around 37 days. The Hipparcos data of BD -08°2823 
includes photometric measurements for 63 epochs only; this 
does not allow a rotation period measurement. 

A different estimation of the rotation period could be ob- 
tained from the HARPS data through the stellar activity indi- 
cators, as shown in Fig. [3] In this figure are plotted the Lomb- 
Scargle periodograms (Press et al. 1992 ) of five HARPS signals: 
the radial velocities, full widths at half maximum (FWHMs), 
contrasts and bisectors of the CCFs, and of the log R' HK activity 
indexes (TableQ]). The fourth last parameters could show signa- 
tures of the stellar activity, and in particular their modulations 
with the stellar rotation; the radial velocity variations could as 
well show signatures of the stellar activity and rotation, but also 
the Doppler reflex motion due to companion(s). The five peri- 
odograms in Fig. [3] all clearly show a peak around 26 days. We 
interpret it as the signature of modulations of the stellar surface 
due to activity (flares, spots, plages...). This kind of phenomena 
are expected for such an active star. As there are sporadic events 
altering the surface of the star and having limited life times, they 
would imply quasi-periodic variations of the shapes of the spec- 
tral lines, with periods near the stellar rotation period, and un- 
conserved phases. The rotation period of BD -08°2823 seems 
thus to be around 26 days from these data. Most of the signals 
show also a small peak near 13-days, which is the first harmonic 
of the main signal; it is not detected in the radial velocities, how- 
ever. In addition, the five periodograms in Fig. [3] show signal 
between 65 and 95 days. It may be the signature of the typical 
duration of cycles of these phenomena. 

The radial velocity periodogram show two extra peaks that 
are not seen for other signals in Fig. [3] at 5.6 and 238 days. The 
fact that these periodic signals are detected only in radial veloc- 
ities suggests it is caused by Doppler reflex motion due to com- 
panions for BD -08°2823 rather than jitter due to stellar activity. 
Moreover, those two periods do not correspond to any harmonics 
nor aliases of the 26 days rotation period. The amplitude of the 
variations (43 m s _1 peak-to-peak, see upper panel of Fig.Q]) and 
the two orbital periods (near 5.6 and 238 days) imply projected 
masses Mp sin i well below the mass of Jupiter. The radial ve- 
locity variations detected by HARPS would thus originate both 
from stellar activity jitter and planetary companions. 
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Fig. 2. Bisector span as a function of the radial velocity (upper 
panel) and of the residuals of the fit with two planets as shown 
in Figs.Q]and[4](lower panel). On the lower panel, the filled cir- 
cles indicate the measurements secured on the 101 -day interval, 
where we fitted the stellar jitter (§ 14.21 and Fig. [6). The ranges 
have the same extends in x- and y-axes on both panels. 



4. A planetary system around BD -08°2823 
4.1. Fit without stellar activity 

We first fitted the radial velocities using a two-planet Keplerian 
model without mutual interactions. The results are plotted in 
Figs.Q]and|4] The inner planet, BD — 08°2823b, produces radial 
velocity variations with a small semi-amplitude K = 6.5 ms _1 , 
corresponding to a planet with a minimum mass Mp sin i 
I - 1 .4 M e , thus similar to the mass of Uranus. Its orbit has a period 
of 5.60 days and is circular or with low eccentricity. The outer 
planet, BD -08°2823c, yields a larger semi-amplitude, namely 
K = 13.4 ms _I ; this corresponds to a planet with a minimum 



mass Mpsin/ = 0.33 Mj up , slightly above the Saturn mass. 
The orbital period is 237.6 days, and the orbit is non-circular 
(e = 0.19). The derived orbital parameters of the system are 
summarized in Table [3] 

The reduced^- 2 of the Keplerian fit is 3.2, and the standard 
deviation of the residuals is <x _ c = 4.3 ms _1 . This is reduced 
when compared to the 1 1-m s _1 dispersion of the original radial 
velocities, but this is clearly higher than the 1.6-ms _1 typical er- 
ror bars on the individual measurements. The extra dispersion, of 
the order of 4 m s ~ 1 , is mainly due to stellar activity j itter as seen 
above. By fitting the two planets without including the stellar ac- 



Hebrard et al.: Planets around the active star BD -08°2823 



5 



BD-082823b 



HARPS 



30 

25 

20 

15 

10 

5 


20 



RV 




FWHM 




10 100 
period (days) 



1000 



10 



> 

cc 



-10 




-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



BD-082823c 



HARPS 



20 



10 



> 
CC 



-10 




-0.1 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 



Fig. 3. Lomb-Scargle periodograms for BD -08°2823 HARPS 
signals. From top to bottom, the five panels show the peri- 
odograms of the radial velocities, of the FWHMs, contrasts and 
bisectors of the CCFs, and finally of the log R' HK activity indexes. 
Each of these five signals shows a peak in the power around 
26 days, which could be interpreted as the stellar rotation period. 
On the upper panel, the radial velocities periodogram shows two 
extra peaks at 5.6 and 238 days, which are not seen in the four 
other panels. The horizontal dashed lines correspond to the false- 
alarm probability of 1 x 10~ 3 . 



tivity, we assume that the stellar jitter would be averaged out over 
the five years of observations, as activity induces quasi-periodic 
effects that are not exactly duplicated in time with the stellar ro- 
tation. This results in a quite large dispersion around the fit with 
respect to the error bars, whereas the periodic signal of the two 
planets remains coherent over the five years of observations. 

Fig. [5] strengthens the interpretation in term of planets of 
the two signals at 5.6 and 238 days. The upper panel shows 
the periodiogram of the radial velocities, exactly as the upper 
panel of Fig. [3] On the second panel of Fig. [5] is shown the pe- 
riodiogram of the radial velocity residuals, after a fit including 
BD -08°2823c only. The standard deviation of the residuals of 
this fit is cr _ c = 5.9 ms _1 . On this periodogram, the peak at 
238 days of course is no longer visible. The main peak is this 
at 5.6 days. Its false-alarm probability is 5.8 x 10~ 5 . It is clearly 
stronger than the peaks due to activity. We note also the pres- 
ence of two other peaks, at 0.8 and 1.2 day, that are the one-day 
aliases of the 5.6-day signal in the frequency space (1 + 1/5.6). 



Fig. 4. Phase-folded radial velocity curves for BD -08°2823b 
(P = 5.60 d, top) and BD -08°2823c (P = 237.6 d, bottom) after 
removing the effect of the other planet. The HARPS radial ve- 
locity measurements are presented with l-cr error bars, and the 
Keplerian fits are the solid lines. Orbital parameters correspond- 
ing to the fits are reported in Table [3] 



We adopted 5.6 days as the actual period of this signal instead 
of one of the two aliases, because it is more likely according to 
the sampling. In addition, its peak is higher than those of the 
two aliases, and Keplerian fits performed with 0.8 or 1.2 day of 
orbital periods produce higher residuals dispersion, and eccen- 
tric orbits. 

Similarly, the third panel of Fig.|5]shows the periodiogram of 
the radial velocity residuals, after a fit including BD -08°2823b 
only. The peak at 5.6 days is removed as well as its aliases, 
and the 238-day signal remains, clearly above the activity peaks. 
Finally, the bottom panel of Fig. [5] shows the periodiogram of 
the residuals after the 2-planet fit shown in Figs. Q] and |4] Only 
lower peaks are remaining, most of them being caused by activ- 
ity. Their false-alarm probabilities are larger than of 1 x 10 3 . We 
note that in addition to the clear peak at one day (which is due 
to the aliases of all the detected signals), all the four panels in 
Fig-E]show a possible peak near 700 days; this could be the sig- 
nature of a third, outer planet. Fits with a curvature in addition to 
the two planets also suggest a possible long-period signal (with- 
out significant effects on the parameters of BD -08°2823b and 
BD -08°2823c). Such long-period signals are not strong enough 
however to claim any detection. Further observations of this tar- 
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Table 3. Fitted orbits and planetary parameters for the 
BD-08°2823 system, with 1-cr error bars. 



Parameters 



BD-08°2823b 



BD -08°2823c 



P [days] 

e 

oj[°] 
K [m s" 1 ] 

T (periastron) [BJD] 
M p sin i [M Jup ] 
Mp sin i [M e ] 
a [AU] 
V r [kms- 1 ] 
cr o-c [ms -1 ] 
reduced x 1 
N 

span [years] 



5.60 ± 0.02 
0.15+0.15 
30 ± 100 
6.5 ± 1.0 
2454637.7 ± 1.6 
0.045 ± 0.007* 

14.4 ±2.1* 
0.056 ± 0.002* 

53.345 ± 
4.3 
3.2 
83 
5.0 



237.6+ 1.5 
0.19 ±0.09 
-233 ±21 
13.4 ±1.0 
2454193 ± 13 
0.33 ± 0.03* 

104 ± 10* 
0.68 ± 0.02* 
0.001 



p. using M t = 0.74 ± 0.07 M e 



get on a longer time baseline are mandatory to establish or not 
the presence of the hypothetic planet BD -08°2823d. 

The lower panel of Fig.|2]shows the bisector spans as a func- 
tion of the radial velocities residuals after the 2-planet fit. This 
does not show an anti-correlation between those two signals, as 
it could be expected in case of spotted stellar surface due to 
activity (see, e.g., Melo et al. 120071 Boisse et al. 2009). The 
absence of anti-correlation, whereas activity is supposed to be 
here the main part of radial velocity residuals, could be due to 
the fact that the rotation velocity is low, or to activity processes 
more complex than simple cool stellar spots. It could also be 
due to the presence of additional companions (as the hypothetic 
BD -08°2823d), that would induce radial velocity variations un- 
correlated with the bisectors. 

We note also that we did not find any anti-correlation be- 
tween the bisector spans and the radial velocities with either one 
or the other planet removed. Such anti-correlation with just one 
planet removed would argue that the other signal is due to activ- 
ity; this is not the case. 

4.2. Fit with stellar activity 

In the previous section we have shown a Keplerian fit of the two 
planets without attempting any fit of the stellar activity jitter sig- 
nal in the radial velocities. Our knowledge of the activity is poor; 
for example the number, locations and sizes of potential stellar 
spots are far to be controlled. We attempt here a naive, phe- 
nomenological approach, based on the fact that activity shows 
signals in the periodograms, so stellar jitter has a periodical na- 
ture that is linked to the rotation of the star. Unfortunately these 
signals are only quasi-periodical: they do not reproduce them- 
selves periodically in an identical form for a long time, as all 
these phenomena have limited life times. Fitting the stellar jitter 
with sinusoids over the 5-year time span do not provide satisfac- 
tory fits. 

We thus tried such sinusoid fits on a shorter time span. We 
chose a 101-day interval, between BJD-2 400000 = 54547.5 
and 54648.5. This interval has a good sampling (27 measure- 
ments), including two 9-night sequences with almost one mea- 
surement per night (8 and 9 measurements in the two sequences) 
at high accuracy (uncertainties between 1.0 and 1.7 ms -1 ). This 
should improve the coherence of the stellar jitter in term of pe- 
riodicity, as well as allowing the 5.60-day-period planet to be 
well sampled. Using sinusoids of periods near 26 days, we fit- 
ted on this 101 -day interval the activity indicators studied in 



E 30 



S 20 




10 100 
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Fig. 5. Lomb-Scargle periodograms of the HARPS radial veloc- 
ities. The upper panel shows the periodogram computed on the 
initial radial velocities, without any fit removed. The second and 
third panels show the periodograms computed on the residuals 
of the fits including BD -08°2823c or BD -08°2823b only, re- 
spectively. The bottom panel shows the periodogram after the 
subtraction of the 2-planet fit. The horizontal dashed lines corre- 
spond to the false-alarm probability of 1 x 10~ 3 . 



Fig. [3] namely the FWHMs, contrasts and bisectors of the CCFs, 
and the logR' HK indexes. The best solution we obtained was 
for a 26.6-day period. Together with the periodograms shown 
in Fig. [3] this allowed us to determine the rotation period: 
26.6 +1.5 days. We fitted on the same interval the radial veloc- 
ities, using Keplerian models for the two planets, and an extra 
sinusoid to approximate the stellar jitter. The parameters of the 
inner planet were free, as the mass of the outer one. The period, 
eccentricity and orientation of the orbit of BD -08°2823c were 
fixed to the values obtained above, the 101 -day interval being 
too short to constrain this 238-day-period orbit. 

The results of these fits agree with those obtained in § 14.11 
The fits are plotted in Fig. [6] which shows on the upper panel 
a good match of the radial velocities and the model, especially 
on the two high-frequency observation sequences that show the 
variations due to BD -08°2823b. The lower panels show the ac- 
tivity signals. The 26.6-day-period sinusoids provide acceptable 
approximations for their variations, which remain however more 
complex than this simple sinusoid shape. The lower panel shows 
the stellar jitter only, i.e. the sinusoid that is added to the two 
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Fig. 6. Fit of the two planets and the stellar activity on a 101 -day 
interval. The upper plot shows the HARPS radial velocities of 
BD -08°2823 with a fit including two Keplerian curves for the 
planets and a 26.6-day-period sinusoid for the stellar jitter. The 
middle plot show the residuals of this fits. The lower plot shows 
on five panels the measurements of the FWHMs, contrasts and 
bisectors of the CCFs, the log/?[j K and the stellar jitter, respec- 
tively from top to bottom. Sinusoid fits with 26.6-day period are 
overplotted. The stellar jitter in the lower panel is the sinusoid 
that is added to the two Keplerian in the upper plot. The vertical 
dashed-line helps visualize the correlations. Error bars at 1 cr are 
plotted for all measurements in this figure. 



Keplerian in the upper panel. It has an amplitude of +5 ms~'. 
The residuals of the radial velocities fit with two Keplerian and 
a sinusoid is plotted in the middle panel of Fig. [6j Its dispersion 
is reduced down to 2.5 ms 4 , which is an improvement by com- 
parison to the 4.3 ms _1 dispersion obtained on the full dataset 
in § 14.11 without an attempt to fit the stellar jitter (or 4.1 ms" 1 
dispersion if computed only on the 101 -day interval, again with- 
out stellar jitter fit). It remains larger than the error bars, because 



of the imperfection of the sinusoid model we used for the stellar 
jitter (and perhaps also because of possible extra planets). We 
restrained this fit to a single sinusoid for the stellar jitter, with- 
out including extra sinusoids at the periods of the stellar rotation 
harmonics. Indeed, these harmonics are not detected in the radial 
velocities (see § 13. 21 and Fig . [3j ■ 

Despite its simplicity, the simple sinusoids allow correlations 
to be seen between the different parameters on this 101 -day time 
span (Fig. [6] lower panels). The FWHMs and contrasts of the 
CCFs are anti-correlated, implying spectral line deformations 
due to activity that let nearly constant their equivalent widths. 
This could be understood as the surface of the CCF is an indi- 
cator of the stellar metallicity (Santos et al. |2002| l. The veloc- 
ity jitter is barely correlated with log R' nK , as well as with the 
FWHM and the bisector of the CCF. The correlation between 
the FWHM and log R' HK is different from what was seen for ex- 
ample in the case of the spotted star CoRoT-7, for which those 
two values are rather anti-correlated (Queloz et al. 120091 ). Also, 
the apparent correlation between log R' HK and the radial veloc- 
ity jitter is different from the picture seen in the cases of the 
act ive sta rs HD 166435 (Queloz et al. l200Tb , GJ674 (Bonfils et 
al. 120071 or HD 189733 (Boisse et al. 2009). And finally, there 
is an apparent correlation between the bisector and the radial 
velocity jitter on this 101 -day time span, possibly with a small 
phase offset (see Fig.|6]l. This possible correlation is also shown 
on the lower panel of Fig. [2] All these relations drawn a picture 
for BD -08°2823 activity that is quite different from a simple 
scenario where the radial velocity jitter is mainly due to dark 
spots on the stellar surface that modulate the shape of the lines 
as the star is rotating. Other phenomena are likely to occur on 
this star; they may include pulsations, convections, flares, plages, 
hot spots... 

The effect of the stellar activity on the observed radial ve- 
locities of BD -08°2823 is of the order of 4 m s . This is lower 
than the semi-amplitudes measured for the two detected plan- 
ets, but not negligible. The error bars reported in Table|3]where 
derived from x 1 variations and Monte Carlo experiments with 
and without stellar activity modeling, as well as from trials and 
errors with different kinds of sinusoids for the stellar jitter. Our 
poor understanding of the stellar activity of BD -08°2823 put us 
to remain cautious on the obtained error bars. The reality of the 
two planets in this range of mass is however well established, 
their signatures being clearly detected at two periods that show 
no signals linked to stellar activity. 

5. Discussion 

Finding transiting planets in the Hipparcos epoch photometry 
annex does not look promising. Our attempt for a-priori de- 
tections did not succeed, and up to now, only two a-posteriori 
detections were performed within the Hipparcos data, in the 
cases of HD 209458b and HD 189733b that were first revealed 
from ground observations. Three main limitations make it dif- 
ficult. First, the error bars on individual photometric measure- 
ments are of the same order of magnitude than the expected sig- 
nal for transits of giant planets, or even slightly larger. Second, 
the sparse time-coverage allows only a few points to be ob- 
tained in a potential transit. These two limitations make any tran- 
sit identification at the limit of detection. Third, stellar activity 
could produce false positives, which are difficult to identify with 
the sparse time-coverage. The ESA mission Gaia (Perryman et 
al. 120011 ) is awaited as the successor of Hipparcos. Its poten- 
tial for planetary transit discoveries in front of bright stars could 
be considered, as its sensitivity and accuracy would be better 
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than Hipparcos. However, the time coverage won't be better 
than that of Hipparcos, which will prohibit well-resolved light 
curve studies. Stellar activity should thus also be a limitation for 
transit detections with Gaia. Dedicated surveys have proven to 
be more efficient for transit detections. Ground-based programs 
allow Jupiter-size planets to be detected, with improving accu- 
racies that now point toward planets with smaller radii (Bakos 
et al. 120091 ), whereas space-based programs allow planets with 
even smaller radii to be detected (Leger et al. 120091 Queloz 
et al. 120091 1, as well as planets on longer periods (Moutou et 
al. 2009b). If accepted, the ESA-proposed space mission PLATO 
(Catala et al. 2009) could permit the detection of such kind of 
transiting planets in front of brighter stars. 

If this search in the Hipparcos data did not provide any de- 
tection of new transiting planets, it allowed the serendipitous dis- 
covery of a new planetary system around BD -08°2823. This 
target was first identified as presenting promising radial velocity 
variations, in which we identified the signature of two new extra- 
solar planets thanks to an intensive monitoring with HARPS. 
The inner planet, BD -08°2823b, is a hot-Neptune with a min- 
imum mass of 14.4 M ffi . Its 5.60-day orbit could be circular or 
slightly eccentric. The outer planet, BD -08°2823c, is slightly 
more massive than Saturn, with a minimum mass of 0.33 Mj up . 
Its orbit is moderately but significantly eccentric, and has a pe- 
riod of 237.6 days. As the masses are low and the orbits are 
distant and nearly circular, the mutual interactions between the 
two planets are negligible. 

The reflex motions that these two planets induce to their host 
star have semi-amplitudes of 6.5 and 13.4 ms 4 , which can be 
distinguished from the 4-ms jitter due to stellar activity. We 
summarize here the arguments that allow us to conclude that 
these two signals are due to planets and not to stellar activity: 

- the signals with periods of 5.60 and 237.6 days are seen only 
in the radial velocities, and are not seen in the shapes of the 
lines nor in the activity indexes; 

- these two periods do not correspond to the rotation period of 
the star, nor to its harmonics or its aliases; 

- with the available data, the signal is coherent over the 5-year 
time span of the observations. 

The radial velocities can be fitted on a short time span by 
two Keplerian and a sinusoid that models the stellar activity. 
Acknowledging that we do not have a good understanding of 
the stellar activity processes, we can alternatively fit the radial 
velocities using two Keplerian only. By doing that we assume 
that the stellar jitter is damped and averaged out over the 5-year 
time-span, and would only lead to a dispersion of the residuals 
of the fit that are larger than expected from the accuracy of the 
radial velocity measurements. Additional planetary companions 
are possible in this system, so the monitoring of BD -08°2823 
should be continued. 

Most of the Neptune-mass and Super-Earth planets found by 
the HARPS GTO program are detected from a survey of about 
200 non-active solar-type stars. They are monitored with numer- 
ous high-precision measurements, using simultaneous thorium 
calibration and exposures of 15-minute duration to average out 
the stellar oscilla tions (Santos et al. I2004bl Udry et al. [2006, 
Lovis et al. [20061 Bouchy et al. |2009l Mayor et al. |2009l . This 
is not the case for BD -08°2823b, that was discovered from 
lower signal-to-noise exposures, without simultaneous calibra- 
tion, and as a serendipitous result of a program that preferen- 
tially selected active stars. BD -08°2823 is one of the most ac- 
tive stars around which a low-mass planet have been found, to- 
gether with CoRoT-7 (Queloz et al. [2009) and a few M-dwarfs 



(Bonfils et al. 2007 Forveille et al. 2009). Those cases support 
the fact that active stars should not be neglected in exoplanet 
hunts, even when searching for low-mass planets. 
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Fig. 7. Stellar metallicities as a function of the planetary masses, 
for 304 known extra-solar planets (top), and for the 96 ones 
among them with orbital periods shorter than 20 days (bot- 
tom). BD -08°2823b and BD -08°2823c are indicated by the 
filled circles. Giant planets are preferentially hosted by stars with 
higher metallicities than low-mass planets. 



BD -08°2823 has a metallicity that is slightly lower than that 
of the Sun. This strengthens the fact that whereas giant planets 
are more frequently found around stars with supersolar metal- 
licities, this tendency is not seen for Neptune-mass nor Super- 
Earth planets. This is shown in Fig. [7] that presents the minimum 
masses of known exoplanets as a function of the metallicity of 
their hosting stars. The data are taken from the compilation of 
the Extrasolar Planets Encyclopedia^. While planets with masses 
larger than 1 Mj up are more numerous around overmetallic stars, 
there are clearly more numerous around undermetallic stars for 
planetary masses below 20 M e . The same effect is seen if we 
select only the planets with orbital periods of 20 days or shorter, 
i.e. in the period regime where most of the low-mass planets are 
detected up to now (lower panel of Fig. [7]). Indeed, among the 
~ 30 planets with a projected mass lower than 0.1 Mj up , only 
four have orbital periods longer than 20 days (namely G1581d, 
HD 40307d, HD 69830c and d). 
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Planet formation models based on the core accretion scenario 
(Mordasini et al. 2009a 2009b) have shown that the sharply ris- 
ing probability of detecting giant planets with stellar metallic- 
ity can be at least partially accounted for by the fact that metal 
rich systems favor the formation of massive planets and that the 
radial velocity technique is most sensitive to massive bodies. 
Observations also seem to indicate that this correlation vanishes 
for Neptune-mass planets, a trend also found in population syn- 
thesis calculations (Mordasini et al, in preparation). These cal- 
culations even show that the correlation reverses for very small 
planetary masses (< 3 - 5 M ffi ), a prediction that will have to 
be confirmed by observations. In fact, a careful analysis of the 
theoretical models indicate that the critical parameter is the over- 
all mass of heavy elements rather than the metallicity. Since this 
mass is determined by the metallicity and the mass of the initial 
proto-planetary disk (which is unknown), it makes a straight in- 
terpretation more difficult. For example, models predict that low 
mass objects orbiting metal rich stars or relatively massive plan- 
ets orbiting metal poor stars are also possible albeit they should 
be rare. 

The BD -08°2823 system also agrees with the tendency of 
low-mass planets to be preferentially found in multiple systems. 
Almost 70 % of the ~ 30 planets with projected masses lower 
than 0. 1 Mj llp are found in multiplanet systems, whereas this ra- 
tio is only ~ 25 % for all the ~ 350 known planets. Considering 
again only the known planets with orbital periods shorter than 
20 days, less than ~ 20 % of them are found to be in multiplanet 
systems: those ~ 20 % are mainly low-mass planets. 

No photometric search for transits have been managed for 
BD -08°2823 from follow-up ground-based observations; de- 
pending on the unknown inclination i of the orbit, the transit 
probability for BD -08°2823b is about 9 %. Its expected depth 
is out of reach of the Hipparcos photometric accuracy. 
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Appendix A: Selection of targets for transiting 
planets search in the Hipparcos database 

We present in this Appendix the systematic search we managed 
in the Hipparcos epoch photometry for planetary transit candi- 
dates (§|2). 

First, we selected 23 304 stars among the 118 204 of the 
Hipparcos epoch photometry catalog, according to the follow- 
ing criteria: 

- B - V > 0.4 (F2 type and later); 

- n > cr n (defined parallax); 

- empty H48 field (reference flag for photometry); 

- H52 field different from D, P, or R (variability types); 

- R+ < 2 R (small stellar radius). 

Then we kept only targets with at most one epoch brighter 
than 3 cr from the average magnitude (in order to remove re- 
maining variables) and at least 40 different available epochs. We 
performed the periods' research on the 17 800 remaining stars. 

For each target, we took the two faintest epochs and we scan 
all the possible periods in order to have those two points in a 
transit. Then we adopted the period producing the lowest flux in 
the drop, and we quantified this solution with the parameter a, 
defined as: 

<V m > - < V oul > 
a =—F 

^ 2 <v m >+<v 0M > 

where < Vi„ > and < V out > are the averaged magnitudes of 
the epochs respectively in and out the drop, and cr the standard 
deviations. The higher a is, the deeper and more significant the 
drop is. We also computed the y 2 f° r the fit of the epochs in 
the drop with a transit curve for a planet with a radius R p = 
(0.1 1 + 0.04) R Q , implying a AVo.n drop in magnitude: 

2 = y (Vi- < Vout > -AV .n) 2 

£—1 rp + rp + rp 

transit U V< u <V 0Ut > ^ u AV .n 
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We finally computed the goodness of fit (gof) for this^ 2 . The 
selection of candidates is based on those two parameters: a and 
the goodness of fit. 
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Fig.A.l. Threshold tuning simulations. The parameters a and 
a' characterize the depth and the significance of the transit, and 
gof the goodness of fit. These parameters are defined in the text. 
Top: Constant stars. Middle: Constant stars (black) and stars with 
transits (red). Bottom: Actual data (blue) and simulations for 
constant stars (black) and stars with transits (red). 



We performed two simulations in order to tune the thresh- 
old: for the 17 800 stars selected above, we assumed in the first 
simulation that all the stars are constant, whereas in the second 
simulation, we assumed that all the stars host a transiting planet 
with a radius chosen randomly within 0.07 - 0.15 R Q . We kept 
the individual error on the photometry on each epoch. 

The upper panel of Fig. |A. 11 shows the goodness of fit as a 
function of a for the first simulation. The distribution is skew 
as a and the goodness of fit are correlated. In order to make 
easier the selection, we fitted the distribution with a line, and 
computed a', which is corrected (with a first-order polynomial) 
from this skew. The middle panel of Fig. IA. II shows the good- 
ness of fit as a function of a' for the two simulations, a' > 5.3 
appears to be an appropriate threshold for transit candidates se- 
lection: indeed, 127 and 489 stars are above this limit in the 
first and second simulation, respectively. We can thus estimate to 
(489 - 127)/ 17800 = 2% the detection rate. With 0.06% chance 
that a given star harbors a transiting extra-solar planet, the 17800 
targets of our latest sample should include ~ 1 1 transiting plan- 
ets. The expected transiting planet detection number with our 
method should thus be slightly less than unity (1 1 x 2% = 0.2). 

The actual data for the 17800 selected targets are plotted in 
the lower panel of Fig. I A. 1 1 The stars with high a' and high 
goodness of fit are eclipsing binaries that were not identified in 
the Hipparcos catalog. In order to remove those binaries, we re- 
moved targets with goodness of fit larger than 2.5. We also re- 
moved targets with R, > 1.3 R Q , and obtained a list of candi- 
dates, sorted by decreasing a'. 

Targets referenced in SIMBAD as binaries, active, or vari- 
able stars were removed from the obtained list, as well as tar- 
gets with already known planets at the time of the observ- 
ing run (HD 70642, HD 39091, HD 10647, HD 4208, HD 17051, 
HD 13445, HD 75289, and HD 83443 - we checked that the 
Hipparcos photometric data did not include periodic variations 
at their periods). 

Using HARPS, we performed follow-up observations of 194 
of these selected targets in December 2004, as part of the pro- 
gram 074.C-0364. The apparent magnitudes of the observed 
stars range from 4.9 to 11.5. Exposures of typically a few min- 
utes duration were obtained, allowing 70 to 90 targets to be 
observed each night. Errors on the measured radial velocities 
are typically of the order of 2m/s. In order to identify hot 
Jupiters, targets were removed from the candidate list after one, 
two, or three observations, according to the following observa- 
tional strategy: 

- observations were stopped after one HARPS measurement 
if: 1) there are two peacks in the CCF; only one SB2 was 
found during our observations (HD 23919), binaries hav- 
ing been previously removed from the candidate list from 
Hipparcos flags and SIMBAD checks; 2) the full width at 
half maximum of the CCF is larger than 15 km/s, which pro- 
hibit accurate radial velocity measurements (17 stars); 

- observations were stopped after two HARPS measurements 
if: 1) the radial velocity difference ARV between the two 
measurements is larger than 2.5 km/s (no SB1 were found); 
2) ARV <5m/s (113 stars); 

- observations were stopped after three HARPS measurements 
if ARV <20m/s(26 stars). 

Radial velocity variations larger than 20 ms~' were mea- 
sured for 37 stars, i.e. 19% of our observed sample. However, 
those variations are mainly caused by stellar activity (see §|5). 



